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Abstract 

A random Eulerian triangulation is a random triangulation where an even number of 
triangles meet at any given vertex. We argue that the central charge increases by one 
if the fully packed 0(n) model is defined on a random Eulerian triangulation instead 
of an ordinary random triangulation. Considering the case n — > 0, this implies that the 
system of random Eulerian triangulations equipped with Hamiltonian cycles describes 
a c = — 1 matter field coupled to 2D quantum gravity as opposed to the system of usual 
random triangulations equipped with Hamiltonian cycles which has c = —2. Hence, 
in this case one should see a change in the entropy exponent from the value 7 = — 1 



to the irrational value 7 = 1 g 13 = —0.76759... when going from a usual random 
triangulation to an Eulerian one. A direct enumeration of configurations confirms this 
change in 7. 
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1 Introduction 



In order to describe the thermodynamic properties of geometrical objects like polymers 
and membranes, it is important to be able to enumerate the possible configurations of 
such objects. One class of problems concerns the folding of these objects. The simplest 
non trivial examples are the compact folding of a self-avoiding polymer, in connection 
with protein folding, and the folding onto itself of a phantom (non self-avoiding) poly- 
merised membrane. Traditionally one has studied the folding statistics of self-avoiding 
polymers by embedding them on a regular two-dimensional lattice and identifying the 
possible folded states of the polymer with the Hamiltonian cycles on the lattice 0. A 
Hamiltonian cycle is a closed curve which visits each vertex of the lattice once and only 
once. Similarly, to describe the possible folded states of a phantom polymerised mem- 
brane one can approximate the membrane by a regular two-dimensional lattice and 
identify the possible folded states of the membrane with those of the lattice. In order 
to facilitate our discussion, we shall restrict ourselves to considering three-coordinate 
or equivalently triangular lattices. The problem of counting the number of folded states 
of the regular, triangulated lattice has been proven to be equivalent to a certain three- 
colouring problem, namely the problem of colouring the links of the lattice with three 
different colours so that no two links which belong to the same triangle carry the same 
colour . This three-colouring problem is a classical mathematical problem which was 
solved by Baxter in 1970 ||. It can also be described as the problem of solving the 
fully packed 0(n) model for n = 2 ||. Similarly, the problem of counting the number 
of Hamiltonian cycles on a regular three-coordinate lattice is equivalent to solving the 
fully packed 0(n) model in the limit n — > 0. This model is critical and describes a 
conformal field theory with central charge c = — 1 [f|, ||. 

However, the use of a regular lattice to describe folding problems is clearly limita- 
tive. For instance, most membranes in nature are fluid rather than polymerised, and 
their modelling requires a random lattice instead. Similarly, one might speculate if the 
complex dynamics of polymers does not call for a random lattice rather than a regular 
one. Making use of the so-called fully packed 0(n) models on a random lattice (cf. 
section |3|) it is possible to generalise the above mentioned folding problems to a random 
lattice. There are, however, some subtleties involved in this process. This is because 
the full packing constraint is very sensitive to the local geometry of the lattice. For 
instance when one solves the Hamiltonian cycle problem on a random three-coordinate 
lattice one finds that c = —2 || |7|, |§. We will explain the reason for this difference 
between the regular and the random case and argue that if we restrict the class of 
random triangulations considered to so-called random Eulerian triangulations we will 
find again c = — 1. By random Eulerian triangulations we mean random triangulations 
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where each vertex is shared by an even number of triangles. Such triangulations also 
appear when we try to generalise the membrane folding problem to a random lattice. 
If we consider simply the fully packed 0(2) model on a random lattice we do not get 
a model describing the folding of a random triangulation. We do also not get a model 
that describes the generalisation of Baxter's three-colouring problem to a random lat- 
tice. In order to have a model which describes the generalised three-colour problem one 
must consider a fully packed 0(2) model where the length of the loops (appearing in 
the graphical representation of the model, cf. section |J) is restricted to being even ||. 
Still, on a random lattice the edge-three-colouring problem is not equivalent to the 
folding problem. In order to have a model which describes the folding problem one 
must consider the edge colouring problem on a random Eulerian triangulation || 
or equivalently the fully packed 0(2) model on a random Eulerian triangulation (note 
that the condition of the loops of the 0(n) model having even length is automatically 
satisfied on an Eulerian triangulation). Eulerian triangulations can also be described as 
triangulations permitting a three-colouring of their vertices |10]. The folding problem 
on a random lattice can hence be viewed as a double three-colouring problem. 

In section |2] and |] we review the properties of the fully packed 0(n) model on 
a regular and on a random lattice respectively and argue that on a random lattice 
we should see a change in central charge for the fully packed model if we change 
the class of triangulations considered from ordinary random triangulations to random 
Eulerian triangulations. In section [| we specialise to the case n — ► and show that 
Hamiltonian cycles on a random Eulerian triangulation should provide us with an 
explicit realisation of a random surface model with an irrational value of the entropy 
exponent, 7. Section |5| describes a numerical solution of the Hamiltonian cycle problem 
on a random Eulerian triangulation and section |6] contains the results - results which 
support our arguments. Finally, in section [7| we conclude and comment on the similarity 
between our Hamiltonian cycle problem and another combinatorial problem which has 
recently attracted a lot of attention, namely the Meander problem 



2 Dense versus full packing of the O(n)-model on 
a regular lattice 

Let us first consider a regular three-coordinate lattice, the honeycomb lattice, and let 
us associate to each vertex, i, of the lattice a n-component spin vector Si with Sf = n. 
The 0(n)-model partition function is then given by 
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where Il<fc,z> is the product over nearest neighbour vertices (TT2] . Expanding the prod- 
uct and integrating over spin variables it can be seen that the only terms which survive 
are those to which one can associate a collection of closed, self-avoiding and non- 
intersecting loops living on the lattice and that the partition function can also be 
written as []I2]1 

SCO = E 0=0 r^ }) =£% } . (2.2) 

{£} Vi 7 {£} 

Here the sum is over all loop configurations {£} having the above mentioned properties, 
V({£}) is the total number of vertices visited by loops and J\f({C}) is the number of 
loops in {£}. The representation (|2.2|) makes it possible to extend the definition of 
the 0(n)-model to non- integer and to negative values of n. It is well-known that for 
n G [—2,2] the 0(n)-model has a second order phase transition at some critical point 
T = T c (n) [[Tj|. The model is also critical for T < T c (n) |H| and in this region of 
the coupling constant space it is denoted as the densely packed loop model, the name 
referring to the fact that the average length of the loops diverges while the number 
of vertices not visited by a loop stays finite (but non zero). The densely packed loop 
model describes a conformal field theory with central charge, Cd, related to n in the 
following way [ 14}| 

6z/ 2 

c d = l-- — -, n = 2cos(z/tt), ue[0,l]. (2.3) 



If we set T = in ( |2.1|) or (|2.2j ) we obtain the fully packed loop model. In this case 
only configurations where all vertices are visited by a loop contribute to the partition 
function. The fully packed loop model has been shown to describe a conformal field 
theory with central charge, Cf given by j|, || 

Cf = c d + l = 2 , n = 2cos(^7r), v E [0, 1 . (2.4) 

1 — v 

If we send n to zero in the fully packed loop model only configurations with one single 
loop visiting all vertices exactly once survive. In this case our partition function thus 
counts the number of Hamiltonian cycles on the honeycomb lattice. According to (|2.4|) 
this statistical mechanical model has central charge c/ = —1. This n — > limit 
corresponds in practice to keeping the linear term in n of the partition function ( |2.2|) . 
For n identically equal to zero, one has c/ = Cd = 0. 

An explanation of the difference in central charge between the densely packed and 
the fully packed 0(n) model on the honeycomb lattice was given by Blote and Nien- 
huis fU]. The fully packed loop model can be viewed as a superposition of a densely 
packed loop model and a SOS-model with central charge c = 1. The extra SOS degree of 
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freedom is a height variable which emerges only in the case of full packing. This height 
variable, h v , lives on the vertices v of the dual triangular lattice. It takes integer values 
and is defined by the demand that for nearest neighbour vertices v and v', \h v — h v / \ = 1 
if the link between v and v' has its dual link occupied by a loop and | h v — h v > \ = 2 
if not [0 . These local rules fix the value of the height variable without ambiguity on 
the entire lattice, up to a global additive constant and a global reversal of all the signs 
of (h v — h v >). We note that the assignment of height variables divides the vertices of 
the triangular lattice into three different sets in such a way that no two neighbouring 
vertices belong to the same set: if we denote the three sets as Si, i = 1, 2, 3 the vertices 
belonging to Si can be characterised by having h v = i (mod 3). Therefore, the possibil- 
ity of constructing the SOS height variable without frustrations essentially relies on the 
fact that the triangular lattice is vertex-three-colourable, i.e. can be divided into three 
sub-lattices of different colour with any two adjacent sites in different sub-lattices. The 
SOS model defined here is equivalent to the zero-temperature anti-ferromagnetic Ising 
model on the triangular lattice |15| and has central charge c = 1 |16[ . A nice way to 
visualise the height variables h v is to mark only those links of the triangular lattice 
which have \h v — h v >\ = 1. This leads to a picture of a 3D-piling of cubes, whose surface 
profile is precisely described by the heights h v , see figure [p. 




h+2 



h+1 




o 



□ 



Figure 1: Construction of the SOS degree of freedom from the fully packed loops, on 
a regular triangular lattice. The thick lines indicate links whose dual link is occupied 
by a loop. The height variable h, defined on the vertices of the lattice, is such that the 
difference of height is ±1 for nearest neighbours connected by a thick line and ±2 for 
nearest neighbours connected by a dashed line. The sign is fixed by demanding that 
h and the number i of the corresponding sub-lattice Si are equal modulo 3. Viewing 
the thick line configurations as a 3D-piling of elementary cubes, the height variable 
corresponds precisely to the height of the surface of the piling. 
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3 Dense versus full packing of the O(n)-model on 
a random lattice 



By duality we can view the 0(n)-model on the honeycomb lattice as being defined 
on a regular triangulation and with that interpretation the model can be coupled to 
quantum gravity by the standard recipe. This leads to the so-called 0(n)-model on a 
random lattice whose partition function reads W7 



Z(A,T) = Ee-^Ei%. (3-1) 
t~S 2 {£} °niA.U 

Here we sum over all triangulations, r with the topology of the sphere S 2 and over 
all loop configurations {£} where it is understood that the loops live on the lattice 
dual to the triangulation. The quantity C T ({£}) is the order of the automorphism 
group of the triangulation r with the loop configuration {£}, N T is the number of 
triangles in the triangulation and A is the cosmological constant. Let us denote the 
triangles traversed by loops as decorated triangles and those not traversed by loops as 
non-decorated triangles. Then we can also write the partition function as 

Z(X, k) = £ e"™- £ -^-^ k»* n^», k = e"* ± (3.2) 

where N n d is the number of non-decorated triangles and Nd is the number of decorated 
triangles. In the coupling constant space (A, k) of the model ( |3.2j ) there is a line of 
critical points beyond which the partition function does not exist. On this critical 
line there is a particular point (A*,/c*) where a phase transition takes place. This 
phase transition is the analogue of the phase transition at T = T c (n) seen on a regular 
lattice. For k < k* the singular behaviour of the partition function is due to the 
radius of convergence in A being reached while for k > k* the singular behaviour of 
the partition function is due to the radius of convergence in k being reached. If we 
approach the critical line from the region where k > k* we reach the densely packed 
loop model on a random lattice. This model has a scaling behaviour characteristic of a 
conformal matter field with central charge c = q coupled to quantum gravity [[H], [18|] . 

If we set e~ x equal to zero while keeping k finite we obtain the fully packed loop 
model on a random lattice. This model has the same scaling behaviour as the densely 



packed model on a random lattice [0. In particular, if we send n — > in the 
fully packed model we get a model describing Hamiltonian cycles on a random three- 
coordinate lattice and this model has c = — 2 || [7]. || . 

The fact that one does not see any change in the central charge when going from 
the densely packed to the fully packed 0{n) model on a random lattice is not in 
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contradiction with the above mentioned argument of Blote and Nienhuis. As we already 
noticed, in order for the construction of the extra SOS degree of freedom to work it 
is necessary that the vertices of the triangular lattice can be divided into three sets 
so that no two neighbouring vertices belong to the same set. A necessary condition 
for this requirement to be fulfilled is that any vertex of the triangular lattice is shared 
by an even number of triangles. On a random lattice a vertex can be shared by any 
number of triangles so the condition is not met, thus preventing the construction of 
the SOS height variable. 

The above requirement of vertex-three-colourability defines a particular class of 
triangulations which we shall call Eulerian triangulations. Such triangulations have 



been studied in ||10|| , in connection with the problem of folding of random lattices. 
For triangulations with spherical topology, the three-colourability requirement and 
the constraint that any vertex is shared by an even number of triangles are actually 
equivalent |19 |. The latter condition can be rephrased by demanding that the number 
of edges leaving any vertex is even, which is the well known property ensuring that the 
triangulation can be drawn in one path without lifting the pen, i.e. can be equipped 
with a Eulerian cycle passing each link once. An example of Eulerian triangulation is 
shown in figure 0. 







£3- 




Figure 2: An Eulerian triangulation made of 8 triangles. Each vertex is shared by an 
even number of triangles. We have indicated the three subsets of different colours, as 
well as the height variables for a particular choice of Hamiltonian cycle (solid line). 



If the argument of Blote and Nienhuis is indeed correct we should be able to recover 
the shift by one in the central charge by restricting the fully packed model to Eulerian 
triangulations. Indeed, for an Eulerian triangulation of spherical topology equipped 
with a fully packed loop configuration, the construction of the SOS variable can be 
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performed globally without frustration. In the case n — > there is a particularly simple 
way of implementing this restriction to Eulerian triangulations, as discussed below. 

4 Hamiltonian cycles on a random Eulerian trian- 
gulation 

Let us consider the fully packed 0(n) model on a random Eulerian triangulation in the 
limit n — > 0. Its partition function reads (cf. equation (|3.2| )) 

= £ ( E TT^V^E^v^ (4-1) 

where the second sum is over random Eulerian triangulations of spherical topology, 
equipped with a Hamiltonian cycle, TC, and consisting of 2N trianglesB. Viewing these 
objects as two-dimensional quantum space-times decorated by configurations of some 
matter field, it follows from the work of David, Distler and Kawai |2(J that behaves 
as 



Zjv-e^A^- 3 + (4.2) 



where 7 is given by 



c-l-,/(25-c)(l-c) 
7 = V \ 2 -, (4-3) 

with c being the central charge of the matter field. In the case of ordinary random 
triangulations equipped with Hamiltonian cycles one finds! that the exponent 7 takes 
the value 7 = — 1 || 0, |J which is characteristic of a matter field with c = —2. What 
we would like to show is that in the Eulerian case 7 takes the value characteristic of a 
matter field with c = — 1, namely the irrational value 

7(C =-1) = Z1^1. (4.4) 

In order to do that we need to determine in ( |4.1| ). Therefore, let us consider a 
random Eulerian triangulation of spherical topology equipped with a Hamiltonian cycle 
and consisting of 2N triangles. Since the triangulation is Eulerian, we can introduce 
a three-colouring of its vertices and this three-colouring of the vertices assigns to each 
triangle one of two possible orientations, + and — . A triangle is said to have orientation 
+ if one encounters the colours of its vertices in the cyclic order 1-2-3 when moving 
counterclockwise along its edges. If the colours are encountered in the cyclic order 
1-3-2 the triangle is said to have orientation — . Obviously, neighbouring triangles have 



4 Note that not all Eulerian triangulations can be equipped with a Hamiltonian cycle. 
5 In this case the formula (4.2) has logarithmic corrections. 
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Figure 3: An arch configuration representing a planar Eulerian triangulation consisting 
of 18 triangles and being equipped with a Hamiltonian cycle. The + and — signs 
alternate along the straight line. Each arch connects a + to a — . 



opposite orientations and any Hamiltonian cycle drawn on the Eulerian triangulation 
visits alternatively + and — triangles. 

We can now use a particular representation which facilitates our analysis. Let us 
represent the Hamiltonian cycle by a straight line having 2N vertices of alternating 
signs. The vertices represent the triangles visited by the loop and the signs the ori- 
entation of these triangles. As mentioned above, neighbouring triangles necessarily 
have opposite orientation. So far each triangle has only been given two neighbouring 
triangles. To completely specify a triangulation we must associate to each triangle one 
more neighbour. This can be done by connecting the vertices pairwise by arches. To 
ensure the planarity of the triangulation, vertices can only be connected either above 
or below the straight line representing the Hamiltonian cycle and different arches can- 
not intersect. Furthermore, an arch must always connect a + to a — . This results in 
arch configurations of the type shown in figure |3|. Counting such arch configurations 
is equivalent to performing the second sum in equation ( |4.1| ). The symmetry factor is 
automatically taken into account. However, two arch configurations differing only by a 
cyclic permutation of vertices or by a reflection in the straight line represent the same 
triangulation. Hence, if we denote the total number of arch configurations involving 
2N vertices as M§ we have 

M§ = 2 x 2N x Z N ~ e »EN NlE -2 L + q f_\ | _ 

We note that the problem of counting Hamiltonian cycles on an ordinary random 
lattice can likewise be reduced to a problem of counting arch configurations. The arch 
configurations appearing in that problem, however, are simpler than the present ones, 
having no signs associated with the vertices. The number of such configurations, Mat, 
can easily be counted using the Catalan numbers, {cn}- One finds [[?[ |J 

M N = c N c N+1 , (4.6) 
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where 

cn = ( ™ ) ~ 4^ iV~ 3 / 2 as iV^oo, (4.7) 

leading to 7 — 2 = 2 x (—3/2), hence 7 = — 1 as stated previously. In the present case 
where we expect to find an irrational value of 7 it is unlikely that a simple counting 
procedure can be found. In the next section we shall present some analytical consider- 
ations on the counting problem and describe a numerical solution. 



5 Numerical solution 

Our starting point is a vertex configuration of the type described above, i.e. a straight 
line with 2N vertices of alternating signs. We want to determine M§, i.e. the number 
of ways of connecting pluses to minuses using two systems of non-intersecting arches, 
one drawn above the line and one drawn below the line. First, let us note that we have 
the obvious inequality M§ < M N from which we get the upper bound \ie < log(16) (cf. 
equations (|4.6|) and Q4.7|) ). Next, let us derive a lower bound on \ie- For that purpose 
we consider the situation where all arches are drawn, say, above the line. In this case 
the requirement of non-intersection of the arches automatically leads to a configuration 
where pluses are connected to minuses since any arch must enclose an even number of 
vertices and hence connects a vertex at an odd position to a vertex at an even position. 
The number of such one-side non-intersecting arch configurations linking 2N points is 
exactly the iV'th Catalan number, cn- From a given one-side arch configuration we can 
construct a two-side arch configuration by flipping a number of arches to the other side 
of the line. The resulting arch configuration obviously contributes to M§. Since there 
are 2 N ways of flipping N arches, we immediately deduce that M§ > 2 n cn from which 
we obtain the lower bound \xe > log(8). Unfortunately, the flipping construction is 
not exhaustive since it produces only configurations where the upper and lower arches 
do not overlap, i.e. configurations where any vertex below a given, say, upper arch is 
connected to another vertex below the same arch even if the connection is made with 
a lower arch. 

We shall describe below a systematic way of counting all the allowed configurations. 
First, let us note that M§ can be written as 

TV 

M§ = J2°N,k (5.1) 

fc=0 

where CW,ft is the number of two-side arch configurations involving 2N vertices and 
having k, say, upper arches. We have the obvious symmetry 

Cn,u = CW,jv-jfc- (5.2) 
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It is possible to calculate Cn,h explicitly for small values of k, as illustrated in Appendix 
A. For the first values of k, Cn,u read explicitly: 



Cn,o — Cat, 
C NA = Nc N , 




(2N-2)\ 



(3N 3 - 2N 2 + 4JV) , 



(5.3) 




(N-2)\(N + 2)\ 
(2N-3)\ 



—N 6 - 15N 5 + —N 4 - 117iV 3 + —N 2 - 192N 

o o o 



(N-2)\(N + 3)\ 
-N ■ 4 N - 2 . 



Besides the trivial case k — 0, the case k = 1 can be easily understood by noticing 
that a two-side arch configuration with only one arch on top cannot have overlapping 
upper and lower arches. Therefore these arch configurations are exactly the arch con- 
figurations produced by starting from one of the Cn,q = Cn arch configurations without 
arches on top and flipping one of its iV arches to the top. This simple procedure breaks 
down as soon as k > 1 and one has to recourse to a more involved strategy described 
in Appendix A. As k increases, however, one very soon runs into rather severe com- 
plications which limit the computation of explicit formulas to the very first values of 
k. In any case, for any finite k, one obtains the large N scaling Cjv.fc ~ 4^, far be- 
low the lower bound on /i E and we expect that the main contribution to M§ comes 
from terms C Njk with k xs N/2. Unfortunately, it does not seem possible to calculate 
such terms analytically. We are not discouraged by this fact, however, because, as 
mentioned earlier, an irrational value of 7 is unlikely to be reproduced by a simple 
counting argument. 

We shall now describe how one can calculate numerically Cjv,fe for (in principle) 
any N and k using a recursive strategy. Again, let us start from a vertex configuration 
having 2N vertices of alternating signs and let us choose k pluses and k minuses which 
are to be connected above the straight line. The remaining N — k pluses and N — k 
minuses are then to be connected below the straight line. Instead of our original 
diagram we now have two sub-diagrams which are to be equipped by one-side arch 
configurations connecting pluses to minuses. Let us denote the number of ways of 
completing with arches a sub-diagram consisting of k pluses and k minuses as c cri(72 ... (J2fe 
where <7j is +1 or — 1 according to whether the vertex at position % in the sub-diagram 
is a plus- or a minus-vertex. For <7j = (—1)*, c CT10 - 2 ... (T2fc is nothing but the /c'th Catalan 
number, c&. Obviously, c (T1(J2 ... (J2fe is non zero if and only if the sequence of oVs is 
globally neutral. There is no closed formula giving c (71(J2 ... (J2fe for an arbitrary sequence 
of (Tj's. This is to be contrasted with the reverse problem, i.e. counting the number of 
+/— sequences compatible with a given arch configuration, which is simply 2^ for N 
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arches since each arch must connect a plus at its left extremity to a minus at its right 
extremity or conversely. This implies in particular the following sum rule 

EW 2 ... CT2 ,= E 2 N = 2 N c N . (5.4) 

{a} arch configurations 

Still, the c r7icr2 ,,, (72k fulfill a certain recursion relation which can easily be implemented 
in a computer program. To derive this recursion relation, let us consider a typical 
sub-diagram, such as the one depicted in figure |j. In this diagram the first plus can 




+ + + + + 



Figure 4: An example of the decomposition of a sub-diagram. 

be connected to any minus at an even position (if the plus is connected to a minus at 
an odd position it is impossible to complete the diagram without crossings). Placing 
the first arch divides the remaining vertices into two groups, those inside the arch and 
those to the right of the arch. These two groups of vertices can be considered separately 
since there can be no arches connecting them. This means that the c (J1(T2 ... (J2JV fulfill the 
following recursion relation 

N 

(5.5) 

where we have used the convention that cq — 1. This recursion relation generalises the 
well known quadratic relation for Catalan numbers. To calculate M§ numerically we 
now use the following obvious strategy 

• Let k run from to , and for each k, select in all possible ways k pluses and 
k minuses to be connected above. 

• Using the recursion relation compute the numbers of ways of completing the 
upper and lower sub-diagram that arise, and multiply these two numbers. 

• Sum the results for all choices of sub-diagrams and all values of k. 

• Multiply by 2. If iV is even, subtract the k = y-term once. 
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In practice, to make use of the recursion formula ( |5.5|) we convert the string a\Oi . . . a^n 
into a number M via the formula: 

2N , _|_ 1 \ 

A/" = £ (^) • 2*-\ (5.6) 

and then save the number of ways of connecting the sub-diagram in a array under 
position M . Our program pre-computes the number of connections for sub-diagrams 
up till M = 2 20 in order to save computer-time. Then formula ( |5.5|) is only used, if a 
sub-diagram has a higher M than this, and such a sub-diagram is quickly reduced to 
diagrams having M less than 2 20 , where the pre-computed values can be used. 



6 Results 

In table [I] we list M§ as a function of N for N = 1, 2, . . . , 20 as determined by the 
computer algorithm described in the previous section. For comparison we list also 
the corresponding number, Mjv, counting arch systems describing ordinary random 



triangulations equipped with Hamiltonian cycles (cf. equation (O)). Furthermore, in 
table § we list Cjf,k f° r N < 12 (cf. equation ( |5.1|) ). We note that in accordance with 
our expectations (cf. section |5|) we observe that the main contribution to M§ comes 
from the Cjv^'s with k « ~. From the data for M§ we shall seek to extract the 



quantities \xe and 7^ (cf. equation (|4.5|) ) using a variant of a ratio method introduced 
in the study of random surface models in pj (for a general discussion, see |]22|j). For 
a generic M N ~ e^ N N^~ 2 {l + 0(1/ N)}, we have that 



$ - ^(^)=» + °(i;)> (6-1) 



(0) 





From these series of estimates we can get improved series of estimates, for instance the 
series 7)^ given by 

■1)^7^ = 7 + 0(^1), (6-3) 

and similarly for /i. In figures ^| and |^, we plot and 7^^ as a function of N for 
p = 2, 3, 4, 5 for the Hamiltonian cycles on the ordinary random triangulations (dashed 
lines) as well as for the Hamiltonian cycles on the random Eulerian triangulations (full 
lines). The series p$ and 7^ are expected to converge faster the larger the value of 
p. However, the larger p the smaller the amount of data. 
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N 


M% 


Mat 


1 


2 


2 


2 


8 


10 


3 


40 


70 


4 


228 


588 


5 


1424 


5544 


6 


9520 


56628 


7 


67064 


613470 


8 


492292 


6952660 


9 


3735112 


81662152 


10 


29114128 


987369656 


11 


232077344 


12228193432 


12 


1885195276 


154532114800 


13 


15562235264 


1986841476000 


14 


130263211680 


25928281261800 


15 


1103650297320 


342787130211150 


16 


9450760284100 


4583937702039300 


17 


81696139565864 


61923368957373000 


18 


712188311673280 


844113292629453000 


19 


6255662512111248 


11600528392993339800 


20 


55324571848957688 


160599522947154548400 



Table 1: The numbers M§ (random Eulerian triangulations) and Mjv (ordinary random 
triangulations) for N = 1, 2, . . . , 20. 



k\N 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 





1 


2 


5 


14 


42 


132 


429 


1430 


4862 


16796 


58786 


208012 


1 


1 


4 


15 


56 


210 


792 


3003 


11440 


43758 


167960 


646646 


2496144 


2 




2 


15 


88 


460 


2250 


10549 


48048 


214344 


941460 


4085950 


17566032 


3 






5 


56 


460 


3172 


19551 


111584 


602514 


3121020 


15655970 


76559920 


4 








14 


210 


2250 


19551 


147288 


1002078 


6320460 


37614016 


213817902 


5 










42 


792 


10549 


111584 


1002078 


7978736 


57977304 


392238792 


6 












132 


3003 


48048 


602514 


6320460 


57977304 


479421672 


7 














429 


11440 


214344 


3121020 


37614016 


392238792 


8 
















1430 


43578 


941460 


15655970 


213817902 


9 


















4862 


167960 


4085950 


76559920 


10 




















16796 


646646 


17566032 


11 






















58786 


2496144 


12 
























208012 



Table 2: The numbers CV,fc for N < 12. 
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Fig function of iV for p = 2, 3, 4, 5 for Hamiltonian cycles on ordinary 

random triangulations (dashed lines) and on random Eulerian triangulations (full lines). 





I 1 


7 


RT (exact) 


log(16) « 2.772589 


-1 


RT (estimated) 


2.772590 ± 10~ 6 


-1.000000 ± 10~ 6 


RET 


2.313 ±0.001 


-0.76 ±0.02 



Table 3: Final estimates for /i and 7 for Hamiltonian cycles on ordinary random trian- 
gulations (RT) and random Eulerian triangulation (RET) along with the exact results 
for the RT-case. 

In table |3] we list our final estimates for /1 and 7 extracted from the series and 7]^ 
for ordinary random triangulations (RT) and random Eulerian triangulations (RET). 
We also give for comparison the exact result for ordinary random triangulations. 

The ratio method reproduces neatly the exact values for the ordinary random trian- 
gulations. When we make the restriction to random Eulerian triangulations we clearly 
see a difference in the behaviour of fjf$ and 7^ . In the case of fiff the convergence 
is as rapid as for ordinary random triangulations but the asymptotic value is smaller. 
We note that our estimate /x# ~ 2.313 respects the bounds log(8) < /x# < log(16) 
derived in section [|, as it should. In the case of 7^ we observe an oscillatory be- 
haviour, not present for the usual random triangulations. Since the amplitude of 
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Figure 6: 7^ as a function of N for p = 2, 3, 4, 5 for Hamiltonian cycles on ordinary 
random triangulations (dashed lines) and on random Eulerian triangulations (full lines). 

the oscillations decreases as N increases the data still allows us to extract a reli- 
able estimate for 7. This estimate for 7 agrees well with the predicted irrational value 
7# = (— 1 — vT3)/6 ~ —0.768 and definitely differs from the value 7 = —1 characteristic 
of Hamiltonian cycles on usual random triangulations. 



We have argued that one should see an increase of the central charge by one unit when 
one considers the fully packed 0(n) model on a random Eulerian triangulation instead 
of an ordinary random triangulation. Our argument was based on an explanation of 
Blote and Nienhuis why one sees an increase in central charge when moving from the 
densely packed to the fully packed 0{n) model on a regular lattice H (cf. section H). 
We have tested our prediction by numerical analysis in the case n — > which cor- 
responds to Hamiltonian cycles on random Eulerian triangulations and on ordinary 
random triangulations respectively. Whereas the latter system has c = —2 the first 
one should according to our prediction have c = — 1. Our numerical analysis confirms 
the prediction. In particular this means that we have found an explicit realisation of a 
random surface model with an irrational value of the entropy exponent, 7, not at any 
point evoking analytical continuation. In addition, our results support the explanation 



7 Conclusion 
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of Blote and Nienhuis. 

Let us denote by A^Ar(g) the number of Hamiltonian cycles on a lattice with iV 



sites and coordination number q. A mean field argument of H. Orland et al. [23] gives 
that 

1 q 
u H = lim — logA^Ar(g) = -. (7.1) 

N^oo Jy e 

In the present case we find that the number of Hamiltonian cycles on random Eulerian 
triangulations behaves as [exp(2.313)] T ^ 2 where T is the number of triangles or equiv- 
alently the number of vertices on the dual lattice. The number of random Eulerian 

rp 

triangulations consisting of T triangles, on the other hand, is known to behave as v^8 



T0|| • This gives the following average value of ujh 



( 2.313 \ V 2 

N « ( -g-J = 1-124, (7.2) 

which is very close to the mean field value - ~ 1.104. For the honeycomb lattice one 
has uj H = ^ « 1.140 §. 

It would be interesting to extend our analysis of the 0(n) model on random Eulerian 
triangulations to other values of n. The case n = 2 is particularly interesting since, 
as explained in the introduction, this case describes the folding problem of a random 
triangulation onto itself, modelling a fluid membrane. As likewise explained in the 
introduction, this folding problem can be viewed as an edge-three-colouring problem on 
a random vertex-three-colourable triangulation. The system of vertex-three-colourable 



triangulations has been found to have central charge c = |L0| and so has the edge 
three-colouring problem on a usual random triangulation ||. Amazingly, when coupling 
these two models one should get a model with central charge c = 2. 




Figure 7: An example of a Meander configuration with 4 connected components. 

Finally, we would like to make a connection between our problem and another 
problem of arches, that of Meanders. Originally, the Meander problem consists in enu- 
merating the number of pairs of arch systems connecting 2N points both from above 
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and from below and such that the path resulting from connecting the upper and lower 
arches is made of a single connected component. This problem turns out to be equiv- 
alent to the compact folding onto itself of a closed linear self-avoiding polymer. By 
extension, one also considers the possibility of creating several connected components, 
with a weight q per connected component. Figure ^ shows an example of a Mean- 
der configuration with 4 connected components. The partition function one wants to 
evaluate is therefore 

= £? c(a,6) , (7.3) 

a, 6 

where a and b are arbitrary systems of arches connecting 2N points on one side (say the 
upper side for a and the lower side for b) of a straight line, and c(a, b) is the number 
of connected components obtained by connecting these arches. As for our problem, 
one expects Zn{q) ~ N 1 ^ -2 at large N. Except for the particular cases q = 1, 
q = — 1 and q —>■ oo which can be solved explicitly, only numerical estimates for \i 
and 7 are known. As for our system, the Meander problem is a problem of interacting 
arch systems. In our case, the interaction comes from the fact that the upper and 
lower arches have to connect pluses and minuses for two sub-sequences of the same 
alternating sequence. In the Meander problem, the interaction concerns the number of 
created connected components. It is interesting to note that this interaction can also 
be reformulated in terms of sequences of pluses and minuses to be connected by arches. 
Indeed, as shown in (TTTJ (in a slightly different language), the quantity c/ c(a ' b ) can be 
written as 

q c{a,b) _ u ^(P(a,{*})+P(b,{*})-Q(a,{*})-Q(b,{*}) ^ ^ 



{cr} a and b 
conncctablc 



for q = u + 1/u and where the sum is over all sets {cr} of pluses and minuses which are 
such that both a and b connect only pluses to minuses (we say: a and b connectable). 
Here P(a, {a}) (respectively Q(a, {cr})) denote the number of arches in a connecting a 
plus (respectively a minus) at the left extremity of the arch to a minus (respectively a 
plus) at the right extremity. For u—1 for instance, this leads to 

M2) = E E 1 = EC™, (7-5) 

a,£) {<?} a and b jcr} 
conncctablc 



to be compared with equation (|5.4j ) . We thus see here that one of the main difficulties in 
the Meander problem actually consists in the determination of the coefficients c ai ... U2N 
for arbitrary sequences. In this sense, solving our problem would be a first step in the 
solution to the Meander problem. 

Acknowledgements We thank J. Ambj0rn, B. Durhuus and J. Jurkiewicz for useful 
discussions and O. Golinelli for a critical reading of the manuscript. 
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Appendix A 



In this Appendix, we will explain how to compute the coefficients C^^ for the first 
values of k. Let us illustrate our strategy by calculating CV,i in a way which can be 
generalised to larger values of k. We again start from a vertex configuration involving 
2iV vertices of alternating signs. Now we must choose a plus and a minus which are 
to be connected above the line. There are two different situations that we need to 
consider. Either the plus chosen is situated to the left of the minus chosen or the 
plus chosen is situated to the right of the minus chosen. Corresponding to these two 
possibilities we write with an obvious notation 

Cjv,i = C N>1 + C Nl . (A.l) 

Let us now consider the first situation (cf. figure |Sp. By choosing the plus and the minus 

+ - + -©- + -+O+- + - + - 

I II III 

Figure 8: Choosing one plus and one minus divides the diagram into three 
sub- diagrams. 

which are to be connected above the line we split the vertex configuration into three 
sub-configurations, I, II, and III. These three sub-configurations all contain an equal 
number of pluses and minuses and we shall denote such configurations as neutral. If 
we considered instead the opposite situation only the sub-diagram in the middle would 
be neutral. The left-most sub-diagram would have a plus in excess and the right- 
most sub-diagram would have a minus in excess. Such diagrams will be denoted as 
positively and negatively charged respectively. As we shall see later the charge of a 
sub-diagram plays an important role in the counting process. Now the vertices of the 
three sub-diagrams must be connected by arches drawn below the line. A vertex of sub- 
diagram II cannot be connected to a vertex of sub-diagram I or III. However, nothing 
prevents vertices from sub-diagram I from being connected to vertices of sub-diagram 
III and we can consider the union of these two sub-diagrams as one sub-diagram. 
In general a series of sub-diagrams can be connected if the resulting diagram is one 
with alternating signs. Now, we can immediately write down an expression for C^\. 
Instead of calculating directly Cat,*;, however, it proves convenient to calculate first the 
corresponding generating functional Ck(z) = J2n=o ^N,k z N ■ Denoting the number of 
vertices in sub-diagram I as 2i and the number of vertices in sub-diagram II as 2j we 
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have 



oo N-lN-l-i oo N-l 

= EE E ^^i-/=EE(^-ih^i-/ 

7V=0 i=0 j=0 N=0 j=0 

= zC(z)±(zC(z)), (A.2) 



where 



C(z) = t c N z N = 1 (A.3) 



By similar considerations one finds for C N \{z 

oo N-XN-X-i 

c.a Cat „• iz = zd(z) I 2. 

dz 



CT + W = E E E c i ^-i-i^ = zC{z) (z±C{z) ) , (A.4) 

N=0 i=l j=0 



the difference from the previous case being the lower limit in the second sum (which 
is due to the fact that in this case sub-diagram I is positively charged). Thus, in total 
we get 

C l {z)=zj-(zC{zf), (A.5) 

which means that 

C Ntl = Nc N , (A.6) 

in agreement with equation ( |5.3| ). For k — 1, the above calculation is clearly not 
the simplest way to reach the desired result, but it has the advantage that it can be 
generalised to the Cjv^'s with a finite k > 1 even if, as we shall see, one very soon runs 
into a rather severe complication. Let us illustrate this by calculating h {z) In 
this case choosing the four vertices which are to be connected above the line divides 
the vertex configuration into five sub-diagrams I, II, III, IV and V which are all neutral 
(cf. figure It is not possible to connect a vertex of a sub-diagram with an odd 

+-+-©- +e+-e-+0+- +- 

I II III IV v 



Figure 9: Choosing two pluses and two minuses divides the diagram into five 
sub- diagrams. 

number to a vertex of a sub-diagram with an even number, but all other combinations 
are allowed. For C 2 h_+_ (^) one thus has, cf. figure |] 

oo N-2N-2-i N-2-i-j N-2-i-j-k 

^ + 1 Z )= 2 EE EE E {Cj+lC k C N ^j-k-l 

N=0 i=0 j=0 k=0 1=0 

+C N - 2 -j-l Cj Ci - Cj C k Q C N -2-j-k-l} z N . (A. 7) 



19 



The first term in the sum represents the configurations where we allow (but not enforce) 
connections between region II and region IV and forbid connections between region III 
and region I or V, The second term represents the configurations where we forbid 
connections between regions II and IV, allowing connections between region III and 
the regions I and V. There is an overlap between the configurations which contribute 
to these two terms and to remove this overlap we must subtract exactly the third term, 
corresponding to situations where the regions II, III and IV are not connected to any 
other region. In the present case it is straightforward to remove the overlap but for 
higher values of k the overlapping becomes a severe complication. 

As long as one is able to keep track of the overlaps, however, there are simple 
diagrammatic rules by means of which Ck(z) can be constructed. Let us assume that we 
are aiming at constructing c^ 1U2 '" a2k (z) where 0\0-i ■ ■ ■ (T 2 k is a given sequence of k pluses 
and k minuses. Choosing the k pluses and the k minuses which are to be connected 
above the line divides the original vertex configuration into 2k + 1 sub-configurations. 
To calculate C% ia2 '" a2k (z) we must sum over all possible ways of connecting these sub- 
diagrams and if necessary subtract overlapping configurations. Assume that we have 
chosen one particular way of allowing connections between the 2/c+l sub-diagrams. Our 
particular choice has arranged the sub-diagrams into groups Gi of possibly connected 
diagrams. The contribution to C^ icr2 "~ U2h (z) from this particular way of connecting the 
sub-diagrams can be written as 



z k c, 



<71<T2...(72fc 



\\FcXz). (A.8) 



Here c aiU2 „, a2k counts the number of ways of connecting the k pluses and k minuses 
above the line. If the total number of sub-diagrams in a group G is denoted as m and 
the number of positively charged sub-diagrams is denoted as m + , it is easily shown 
that Fq(z) takes the form! 

i Jm—l 

For instance in the case of ~ 2 + ~ one finds immediately 



C+- + -(z) = 2z 2 C{z) 



UzC{z)) 
dz 



2 + ^ 2 C{zf\^C{z)) - 2z*C(zf ±(zC(z)), 



which is easily seen to coincide with the expression ( |A.7| ). Using the rules ([A. 8 



and (|A.9| ) as well as the explicit expression for C(z) (cf. equation ( |A.3|) ) we have 



6 In any given group of connected sub-diagrams the number of positively charged diagrams must 
be equal to the number of negatively charged diagrams. 
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been able to determine Cn,2 and Cjv,3, as given in equation (|5.3|). It is obvious that 
as long as we consider k finite, Ck(z) can be expressed in terms of a finite number of 
terms involving only the function C(z) and its derivatives and thus C^,k ~ 4^ at large 
N. 
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